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Hyperbolic metamaterial (HMM), a sub -wavelength periodic artificial structure with hyperbolic dispersion, 
can enhance the spontaneous emission of quantum emitters. Here, we demonstrate the large spontaneous 
emission rate enhancement of an organic dye placed in a grating coupled hyperbolic metamaterial 
(GCHMM). A two-dimensional (2D) silver diffraction grating coupled with an Ag/Al203 HMM shows 
18-fold spontaneous emission decay rate enhancement of dye molecules with respect to the same HMM 
without grating. The experimental results are compared with analytical models and numerical simulations, 
which confirm that the observed enhancement of GCHMM is due to the outcoupling of non-radiative 
plasmonic modes as well as strong plasmon-exciton coupling in HMM via diffracting grating. 

In view of cutting-edge scientific advances in the field of nanophotonics, unique multi-scale metamaterials have 
been recognized as versatile engineered structures to mitigate classical limitations of optics and enable numer- 
ous challenging applications. The main core in all metamaterials comprises of fabricating a medium composed 
of unit cells of size far below the excitation wavelength. In contrast to existing optical metamaterials in the visible 
and near-infrared wavelength regions, hyperbolic metamaterials (HMMs) are distinguished as promising class of 
metamaterials owing to their characteristic features like hyperbolic form of the iso -frequency dispersion curve 
and the directional nature of light propagation^"^. The artificial HMMs can be analogous to uniaxial metacrystals 
with an extremely anisotropic dielectric tensor, s = diag [Sxx, Syy, Szz] such that Sxx = Syy and • Sxx < 0, which 
leads to a hyperbolic dispersion exhibiting dielectric properties (s > 0) in one direction and metallic nature in 
another direction (s < 0). This distinctive physical consequence of dispersion relation enables the propagation of 
waves with very high wavevectors, known as high-k waves, showing a broadband singularity in the density of 
photonic states^. These interesting properties of HMMs promise a large variety of substantial futuristic applica- 
tions such as spontaneous emission enhancement^"^", nanoimaging^^"^^ negative refraction^^"^^, super-planckian 
thermal emission^^, biosensing^^ and Dyakonov plasmonics^^. The hyperbolic dispersion of HMMs leads to the 
existence of high photonic density of states, which can be investigated by understanding the modifications in 
spontaneous emission of emitters placed in the vicinity of HMM^. 

In cavity quantum electrodynamics, closed cavities (microcavities) and open cavities (photonics crystals) have 
been employed as significant tools to understand cavity quantum electrodynamics and hence for the study of 
modification of spontaneous emission enhancement by the Purcell effect^°"^^. However, the incapability for 
'microcavities' and 'photonic crystal' systems to obtain high quality factor resonance make them unviable for 
Purcell enhancement. This incompatibility of cavity technologies pushed towards the discovery of novel systems 
whose exhibit broadband Purcell effect. Hyperbolic metamaterials are promising and reliable systems in com- 
parison to conventional cavity approaches due to the existence of a large number of electromagnetic states, which 
lead to a divergence in the photonic density of states allowing a broadband control. In particular, HMMs play a 
vital role to confine the power emitted by emitters in large spatial wavevector channels. This will open new routes 
for innovative photonic device applications from broadband single photon sources to strong coupling of emitters 
and plasmonics. In past years, many research groups have explored about modifications in the spontaneous 
emission rates of organic dyes and quantum dots positioned in the close vicinity of the surface of a metamaterial 
with hyperbolic dispersion^"^. However, the Purcell enhancement of HMM is limited at the far-field due to the 
non-radiative behaviour of plasmonic modes in HMM. Recently, nanopatterned hyperbolic metamaterials have 
been proposed for outcoupling the non-radiative plasmonic modes and to enhance the spontaneous emission 
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Figure 1 | Fabrication of GCHMM (a) Schematic of fabricated Ag diffraction grating coupled Ag/Al203 HMM (GCHMM), which consists of 6 pairs of 
Ag/Al203. (b) SEM image of 2D Ag diffraction grating on top of the PMMA layer with an average period of 500 nm and hole size of 160 nm. 



rates of fluorescent molecules^^'^^. Here we propose another novel 
approach, which is based on a grating coupled HMM (GCHMM) 
configuration to further improve the spontaneous emission rate of 
emitters placed nearby HMM, by exploiting the unique property of a 
hypergrating to outcouple and extract the non- radiative plasmonic 
modes. We have recently demonstrated that a grating coupled HMM 
geometry (hypergrating geometry) can be used to excite both surface 
and bulk plasmon modes from a HMM system^^. In the present work, 
we experimentally and numerically investigate the large spontaneous 
emission rate enhancement of fluorescent molecules using a 
GCHMM configuration. 

Results 

Design and fabrication of GCHMM. The designed grating coupled 
HMM (GCHMM) structure (Fig. 1 (a)) is fabricated by following 
sequential depositions of 12 alternated layers of alumina (AI2O3) and 
silver (Ag) thin films by using conventional electron beam and 
thermal evaporation techniques, respectively. All the thin films are 
grown over micro -glass substrates and AI2O3, Ag pellets are used as 
source materials. Spectroscopic ellipsometry measurements have 
been carried out to estimate the thicknesses of grown thin films. 
The measured thickness of AI2O3 and Ag are 23 nm and 12 nm, 
respectively. Therefore, in the present investigation evaluated fill 
fraction of Ag is 34%. To investigate the influence of designed 



HMM on spontaneous emission enhancement a dye doped 
PMMA layer (thickness around 100 nm) is spin coated over a pre- 
deposited AI2O3 spacer (12 nm) on the HMM structure (see 
Fig. 1(a)). In this process, an organic dye (Coumarin 500) is 
initially dissolved in ethanol (0.3% by wt.) solution and in the next 
step resultant solution is dissolved in PMMA resist. For comparative 
analysis, a reference sample is fabricated by spin coating a layer of 
dye-doped PMMA over a pre-deposited AI2O3 spacer (same 
thickness of 12 nm) over glass substrates. The dielectric constants 
of thin films of Ag and AI2O3 are experimentally determined using 
spectroscopic ellipsometry. The real (red line) and imaginary (blue 
line) permittivity values of Ag as well as dielectric constants (black 
line) of AI2O3 are shown in Fig. 2 (a). One can clearly see that real 
permittivity values of Ag show a well-known decreasing trend 
towards strongly negative values. This is due to the Drude-type 
response of free electrons in metals at higher wavelengths^^. In the 
fabricated HMM structure, the individual metal and dielectric layers 
dimensions are satisfying the criteria of effective medium theory 
(EMT) to achieve homogeneity. Therefore the uniaxial dielectric 
tensor components of an anisotropic medium of HMM are 
evaluated using EMT^^. In the calculation, the optical constants of 
Ag and AI2O3 are taken from experimentally obtained values 
(Fig. 2(a)). The evaluated dielectric permittivity tensor components 
of the fabricated HMM (i.e SuiSx = Sy) < 0 and 8 j_ (e J > 0) imply a 
hyperbolic dispersion above the wavelength of 430 nm, at which 
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Figure 2 | Characterization of GCHMM: (a) Experimentally determined permittivity values of Ag and AI2O3 thin films. The red and blue lines represent 
the real and imaginary permittivity values of Ag, respectively and black line represents the permittivity values of AI2O3. (b) Real parts of effective 
permittivity of Ag/Al203 HMM determined with effective media theory. The Ag/Al203 HMM shows hyperbolic dispersion at }^ > 430 nm. 
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Figure 3 | Reflectance spectra as a function of excitation wavelength: (a) Reflectance spectrum of various samples obtained at incident grazing angle of 
50°. The studied samples are shown in Figs, (b) to (e). The GCHMM sample (iv) shows four prominent reflectance dips corresponding to the excited 
plasmonic modes (both surface and bulk plasmon polaritons) from the geometry. 



transition occurs from elliptical to hyperbolic dispersion as shown in 
Fig. 2(b). 

A grating coupling technique is used to excite plasmonic modes 
associated with HMM. By following the electron-beam lithographic 
technique, a 2D diffraction grating is patterned above the dye mixed 
PMMA layer over the HMM, followed by deposition of an Ag layer of 
thickness 20 nm. The fabricated grating has an average period of 
500 nm and average hole diameter of 160 nm. The total area of 
patterned grating structure over HMM is around 5 mm X 5 mm. 
The SEM image of patterned 2D Ag diffraction grating is shown in 
Fig. 1(b). During the diffraction grating fabrication, Ag film of thick- 
ness 20 nm is directly deposited on the PMMA diffraction grating, 
but not directly on the dye mixed PMMA layer. Note that a thin 
PMMA spacer layer is at the interface between Ag film and dye mixed 
PMMA layer because a thin layer of PMMA resist always left during 
the resist development (etching) process (see Fig. 3 (e)). Hence this 
thin spacer layer can avoid the nonradiative contribution due to 
quenching at the Ag-Dye mixed PMMA interface. According to 
grating coupling principle, the surface plasmon modes can be excited 
when the wavevector of the grating diffraction orders are greater than 
the incident light. Under this condition, diffraction orders are no 
longer propagating waves, but evanescent field, and the enhanced 
wavevector of evanescent field is responsible for the coupling of 
incident light to surface plasmon modes by following the condition 
kspp = nokosinO ± rnkg^ ± nkgy. Being 0 the incident grazing angle, Hq 
is the refractive index of incident medium, ko = 2n/l is the vacuum 
wavevector, m and n are the grating diffraction orders and kg = 2n/A 
is the grating wavevector with A being the grating period. On the 
other hand, by introducing a diffraction grating on top of the HMM, 
it is possible to diffract light and produce a wide range of wave vectors 
into the HMM. Due to the existence of impedance mismatch at the 
various openings, the generated wave vectors can couple through the 
surface modes^^. It should be noted that the large thickness (around 
100 nm) of dye dissolved PMMA layer can decrease the coupling 
efficiency of the proposed structure. Nonetheless, unintentional scat- 
terers are common in dye dissolved PMMA layer and these scatterers 
also enable the excitation of surface plasmon modes through the 
coupling of incident light to the surface plasmon polaritons 
(SPPs)^^'^°. These two phenomena contribute to the excitation of 
high-k modes inside the HMM. Since the emitter is placed inside 
the GCHMM, it is a challenging operation to couple the incident 
radiation to the emitter because of the impedance mismatch between 
air and HMM. However, the grating coupling helps to outcouple the 
highly confined modes from the structure to the far field. In order to 
show the excited plasmonic modes from HMM, the reflectance spec- 
trum of different samples for p-polarization is obtained and shown in 
Fig. 3 (a) (also see Supplementary Fig. 1 & Fig. 2). Here incident 
grazing angle is set to be 50°. One can observe that 2D Ag grating 



coupled HMM supports four narrow reflectance minima that corre- 
spond to expected plasmonic modes (both surface and bulk plasmon 
modes). The reflectance spectrum of GCHMM is simulated using 
scattering matrix method (SMM)^^'^^ that is shown in Fig. S3 
(Supplementary Fig. 3). A good agreement with experimental result 
is obtained. In addition, Ag film deposited HMM shows a broad 
reflectance minimum at higher wavelength (750 nm), which is due 
to the presence of Ag scatterers intermixed at the interface with dye 
doped PMMA layer as experimentally verified (see Supplementary 
Fig. 4). 

Experimental results. Time-resolved photoluminescence measure- 
ments are performed to study the decay rate of various samples. An 
ultrafast optical setup has been used (see Supplementary Fig. 5); it 
consists of a Ti: Sapphire tunable femtosecond laser, a Pulse Picker, a 
Second Harmonic Generator and a spectrofluorometer used for 
time-correlated single photon counting (TCSPC) measurements. 
The time resolution of the TCSPC instrument is <5 ps. The 
Coumarin molecules are excited by using a pulsed laser at 380 nm 
with a pulse width of about 120 fs and a repetition rate of 4 MHz. 
The maximum emission wavelength of Coumarin 500 dye dissolved 
PMMA is observed at 470 nm for 380 nm excitation wavelength (see 
Supplementary Fig. 6). In order to show the transition from elliptical 
to hyperbolic dispersion, the investigated emission wavelengths are 
varied. The short-living excitonic states of the emitters placed in the 
vicinity of the HMM and the measured fluorescent lifetime as 
a function of emission wavelength, would represent a clear 
signature of the transition from elliptical to hyperbolic dispersion. 
The fluorescence time decay curves of four samples at spectral 
regions such as elliptical dispersion (420 nm), the critical 
wavelength (430 nm), and the hyperbolic dispersion (450 nm) are 
shown in Fig. 4. Four different samples such as reference (pink 
curve), HMM (red curve), Ag film deposited HMM (black curve) 
and GCHMM (blue curve) are investigated for comparison. It was 
not possible to fit the decay rate of all four samples using a single 
exponential function because the detected signals arise from the 
collective response of molecules randomly distributed in PMMA^^. 
Therefore the data have been fitted using three exponential 
functions, R(t) = Bie~*^^' + B2e~^^'' + B^e'^^'' with Xi being the 
decay times. Since longer time (13) is attributed to uncoupled dye 
molecules, located above the coupling distance from the HMM, we 
use shorter decay times (xi and X2) to predict the decay rate 
enhancement, since shorter decay times are related to molecules 
strongly coupled with HMM^^. The transition from elliptical to 
hyperbolic dispersion is evident from the obtained curves, showing 
a large variation in time decay for HMM compared to reference 
sample when the emission wavelength is varied from elliptical to 
hyperbolic region (also see Supplementary Fig. 7 & Fig. 8). The 
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Figure 4 | Time-resolved photoluminescence measurements of Coumarin dye on various samples, Ref (pink), HMM (red), HMM with Ag (black). 



and GCHMM (blue)) with emission wavelength: (a) in elliptical region {X^ 
region (X^ = 450 nm). 

coupling of the high-k metamaterial states is responsible for the 
observed behaviour of HMM^. The observed response of HMM 
definitely supports the existence of high k-modes in the fabricated 
Ag/Al203 HMM. However, the observed results are not in good 
agreement with previously reported results^. This could be due to 
the low internal quantum yield of the Coumarin 500 dye in PMMA, 
broadband dye emission spectrum and higher metallic fill fraction. 
The decay rate enhancement is defined as the inverse of decay time. 

To verify the large decay rate enhancement of proposed config- 
uration, the GCHMM sample is compared with HMM and Ag film 
deposited HMM samples. Throughout the wavelength spectral range 
of the HMM, one can see that the GCHMM (blue curve) shows a 
large decay rate variation with respect to HMM (red curve) and Ag 
film deposited HMM (black curve). In addition, the decay rate vari- 
ation increases when the emission wavelength is varied from the 
elliptical to hyperbolic spectral region of the HMM. The maximum 
decay rate is obtained for GCHMM for all wavelengths, which are 1/ 
0.035 ns~\ 1/0.030 ns~\ and 1/0.031 ns"' for 420 nm, 480 nm and 
510 nm, respectively. The corresponding decay rates at the respect- 
ive wavelengths are 1/0.220 ns~\ 1/0.491 ns~\ and 1/0.575 ns"' for 
HMM and 1/0.071 ns"\ 1/0.353 ns"' and 1/0.534 ns"' for Ag film 
deposited HMM. Note that the first decay time (xj is considered 
here for comparison. This large decay rate of GCHMM is attributed 
to the strong coupling of emitters to HMM via metallic diffraction 
grating. Specifically, the metallic grating excites the plasmonic Bloch 



- 420 nm), (b) at critical wavelength (X^. = 430 nm), and (c) in hyperbolic 



modes (high-k modes) as well as surface plasmon polaritons of 
HMM and there is a strong overlap between these modes with the 
quantum emitters, which leads to a broadband enhancement of 
photonic density of states (PDOS). In addition, 2D diffraction grat- 
ing scatters the high-k modes of HMM into well-defined free space 
modes. The sub -wavelength confinement of the emitter inside the 
GCHMM structure also enhances the decay rate. Furthermore, 
improved decay rate is observed for Ag film deposited HMM at 
certain emission wavelengths as compared to HMM. It is attributed 
to the presence of Ag scatterers intermixed at the interface with dye 
doped PMMA layer and nonradiative contribution due to quenching 
at the Ag-dye doped PMMA interface. 

In order to emphasize the improved plasmon- exciton coupling in 
GCHMM, normalized lifetimes of GCHMM and Ag film deposited 
HMM are compared (the fitting curves are shown in Supplementary 
Fig. 9). The first (xj and second (X2) decay times of both samples, 
normalized with respect to HMM as a function of emission wave- 
length, are shown in Fig. 5(a) and Fig. 5(b), respectively. It is evident 
from those figures that same decay rate variation is obtained for both 
decay times. For GCHMM, lifetime values are higher in elliptical 
region and smaller in hyperbolic region. This behaviour represents 
the preferential emission of high-k modes from GCHMM. However, 
the lifetimes of Ag film deposited HMM are randomly varying along 
the spectral range showing that the presence of Ag film does not 
support the preferential emission of high-k modes. The smaller life- 
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Figure 5 | Lifetimes of dye on GCHMM and Ag film deposited HMM normalized with respect to HMM sample: for (a) first decay time (xi) and (b) 
second decay time (X2). The solid lines are for eye guide. Enlarged plot of normalized decay time of GCHMM is shown in the inset of (a) and (b). 
Normalized lifetime is studied as a function of emission wavelength from elliptical to hyperbolic dispersion. In the case of GCHMM, a maximum decay 
rate is observed at 510 nm for both decay times. 
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Figure 6 | Numerical simulation results (a) Calculated decay rate enhancement of HMM and GCHMM with period 500 nm. Enlarged plot of decay rate 
enhancement of HMM is shown in the inset of (a), (b) Tuning of decay rate enhancement of GCHMM with grating period. 



time values of GCHMM in hyperbolic region are due to the hyper- 
bolic dispersion. It should be noted that the lifetime values of 
GCHMM behave non-monotonically as function of the wavelength, 
they decrease by reaching a minimum value at 510 nm, then in the 
same hyperbolic region show a change of slope towards higher life- 
time values. The observed behaviour is a clear signature of the out- 
coupling ability of GCHMM. According to Xi, a maximum of 18 -fold 
decay rate enhancement is obtained for GCHMM at 510 nm emis- 
sion wavelength, whereas 5 -fold decay rate enhancement is obtained 
for T2. In comparison to Ag film deposited HMM, the maximum 
decay rate enhancement for GCHMM is observed at 510 nm, which 
is around 17 times and 4 times higher for Xi and T2, respectively. The 
observed large spontaneous emission rate enhancement of GCHMM 
is mainly due to radiative recombination rate, whereas non-radiative 
recombination rate is negligible in GCHMM. Thus the observed 
behaviour of GCHMM strongly supports the influence of diffraction 
grating on the spontaneous emission rate enhancement. The main 
advantage of GCHMM compared to the other two samples is the 
ability to outcouple the high wavevector plasmonic modes from 
HMM to far- field, by the detection of fast- decaying signals from 
emitters strongly coupled to the multilayer structure^^. In the case 
of HMM, the field inside the multilayer is highly confined and it is 
unable to be detected at far field due to the evanescent nature of field 
at the top surface of HMM. However, the grating coupler helps to 
translate the evanescent field to a propagating field for far- field 
detection. This is evident from the time-resolved photoluminescence 
measurements. 

Simulation results. A numerical simulation has been implemented 
to obtain a complete analysis of the enhanced emission rates 
observed in the experimental section. In order to study the decay 
rate enhancement analytically, we use a semiclassical approach in 
which quantum emitters have been considered as radiating point 
dipoles^l According to this method, the decay rate enhancement 
with respect to vacuum of a point dipole with dipole moment '/i', 
placed at a distance 'cf from a HMM is given by^'^'^, 

GO 

^ r ^ X 3w r ^ . kiidk 

0 

1^ Hl [(1 + rX - (1 - rp)kl] + nl{l + rp)kl | 

Where r] is the intrinsic quantum yield of the emitter, kg is the 
wavevector of the medium where the emitter resides, ku is the 



wavevector parallel to interface of the HMM, kz = ykl — k\ and r^, 

Yp are the Fresnel reflection coefficients for s- and p-polarization, 
respectively. Here the Fresnel reflection coefficients of the HMM 
are obtained with transfer matrix method^^. In our simulation, the 
relative permittivity of Ag is obtained from Drude model, ^Agio^) = 

1 - (^o)j/((o^ + ivco)^ with (Op =11.5 fs'' and v = 0.083 fs'' '\ Note 

that the point dipole (emitter) resides in PMMA and refractive index 
of PMMA is set to be 1.49 for wavelengths in the range of 400- 
800 nm^^. Here the separation distance and intrinsic quantum yield 
of Coumarin dye is taken as 12 nm and 0.98^^, respectively. 

Figure 6 shows the decay rate enhancement of HMM and 
GCHMM. As shown in the inset of Fig. 6 (a), an improved decay 
rate enhancement is observed for HMM in the hyperbolic spectral 
region (above 430 nm). Due to the strong coupling between the 
emitter and the high-k modes in HMM, the decay rate enhancement 
monotonically increases in hyperbolic region. In particular, this 
coupling strength increases for shorter interaction distances (d/1). 
In comparison to HMM, GCHMM shows improved decay rate 
enhancement for all wavelengths. However, for a particular wave- 
length band, GCHMM provides very large enhancement of P with 
respect to HMM. For 500 nm grating period, the maximum P is 
obtained at 497 nm wavelength that is 12 times higher than HMM. 
It shows that the simulation results are in good agreement with 
experiments. However, there is a slight discrepancy in wavelength 
at which maximum emission rate is observed. This can be attributed 
to the non-uniform periodicity (average period = 500 nm) of fabri- 
cated diffraction grating as compared to uniform periodicity 
(500 nm) considered in the simulation. According to the simulation 
model, the maximum value of p in GCHMM depends on the peri- 
odicity of the diffraction grating. Note that the periodicity of the 
grating is introduced in the calculation of Fresnel refiection coeffi- 
cients for s- and p-polarization, i. e, k^ = nokoSinO ± m(27r/AJ ± 
n{2nlKy) with A^^ = A^ = A for grating with square lattice symmetry. 
The decay rate enhancement of GCHMM with different grating 
periods is shown in Fig. 6(b). It shows the possibility of tuning the 
spectral position of the maximum P and broadband enhancement of 
P by varying the grating period. According to Fig. 6(b), the larger 
enhancement of P is observed for smaller period (A = 475 nm) that 
is due to the fact that smaller periods provide better momentum 
match and outcouple most of the high-k modes from HMM^^. 
Since HMM provides broadband decay rate enhancement in hyper- 
bolic region, a single period diffraction grating can outcouple a single 
wavelength as observed in our simulation. However, up on designing 
a proper chirped grating a broadband outcoupling of plasmonic 
modes is possible^^'^^. Therefore the broadband outcoupling of single 
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Figure 7 | FDTD simulation results for HMM and GCHMM. Cross-sectional map of the electric field distribution in the hyperbolic region with emission 
wavelength such as 480 nm, 497 nm, and 510 nm: (a)-(c) for HMM and (d)-(f) for GCHMM. For GCHMM, the period of the grating is set to be 
500 nm. 



photon as well as outcoupling of greater number of high-k modes can 
be achieved by engineering a modified diffraction grating geometry. 

We further performed finite difference time domain (FDTD) 
numerical simulation to show the effect of grating coupling in 
GCHMM. In order to compare the obtained results, the electromag- 
netic field distribution along the HMM and GCHMM is simulated by 
assuming that the quantum emitter (point dipole) with a specified 
emission wavelength resides in PMMA above the multilayer. In 
hyperbolic region, the electric field (EF) cross-sectional map (xz- 
plane) of both HMM and GCHMM as a function of emission 
wavelength is shown in Fig. 7. According to Fig. 7(a) to 7(c), the 
magnitude of EF for HMM increases when the emission wavelength 
is varied from 480 to 510 nm. This is due to the monotonic increase 
of the fluorescence decay rate of HMM in the hyperbolic region. In 
comparison to HMM, GCHMM provides enhanced EF on the sub- 
strate as well as above the grating (in air), which evidently shows the 
outcoupling effect of GCHMM. It is clearly visible that the EF dis- 
tribution on GCHMM extensively varies with emission wavelength. 
The maximum EF distribution for GCHMM along the multilayer 
and the strong coupling between grating, and point dipole is 
observed for 497 nm emission wavelength. This behaviour of 
GCHMM definitely supports the maximum decay rate enhancement 
observed at 497 nm in Fig. 6(a). Note that the enhanced EF distri- 
bution is observed close to the surface of the HMM in which higher 
plasmonic density of states can be outcoupled. In the case of 
GCHMM, the enhanced field distribution observed at 497 nm defi- 
nitely supports a strong overlap between the excited high-k modes 
and quantum emitter, which can provide a broadband enhancement 
of photonic density of states, in addition to an outcoupling effect 
from GCHMM. As mentioned before, the outcoupling strength 
depends on the periodicity and the geometry of the diffraction grat- 
ing. In short, the outcoupling of non-radiative plasmonic modes as 
well as strong plasmon-exciton coupling through the diffraction 



grating is responsible for the observed large spontaneous emission 
decay rate enhancement in GCHMM. 

Discussion 

In summary, we designed and fabricated a grating coupled hyper- 
bolic metamaterial based on Ag/Al203 multilayers to show the out- 
coupling effect and large spontaneous emission rate enhancement of 
fluorescent molecules. In comparison to Ag/Al203 HMM without 
diffraction grating, we obtained 18-fold spontaneous emission rate 
enhancement of dye molecules placed in the close proximity of 
GCHMM. Then we performed simulation based on a semiclassical 
approach to show the decay rate enhancement analytically. We fur- 
ther carried out FDTD numerical simulation to show the enhanced 
field distribution and outcoupling effect in GCHMM. A good agree- 
ment with experiment and simulation is obtained, which confirms 
that the observed enhancement of GCHMM is due to the outcou- 
pling of non-radiative plasmonic modes as well as strong plasmon- 
exciton coupling in HMM via diffracting grating. The proposed 
grating coupled configuration can offer the possibility of engineering 
the plasmonic density of states for Purcell factor enhancement for 
potential applications including biosensing^°, surface-enhanced 
Raman spectroscopy^^ and single photon sources^\ 

Methods 

Sample fabrication and characterization. Hyperbolic metamaterials were realized 
by sequential deposition of AI2O3 and Ag layers on a glass substrate (Micro slides 
from Corning) using E-beam evaporation of AI2O3 pellets (from Kurt J. Lesker 
Company) and thermal evaporation of Ag pellets (from Kurt J. Lesker Company), 
respectively. The deposition rate of both materials are set to be 0.3 A/s. Coumarin 500 
dye (from Exciton) was initially dissolved in ethanol (0.3% by wt. in ethanol solution) 
and then the solution was dissolved in PMMA resist (950PMMA C2 Resist from 
MICROCHEM). As prepared solution was spin coated on the HMM sample at 
5000 rpm, in order to get a thickness of around 100 nm. The thickness was measured 
using a Stylus Profilometer (KLA-Tencor P-6). Two-dimensional diffraction grating 
on top of the dye dissolved PMMA was fabricated using Electron-beam lithography 
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(Tescan Vega). Initially, MMA resist (8.5 MMA EL 1 1 from MICROCHEM) was spin 
coated on the sample at 4000 rpm and baked at 180°C for 2 min. After some time, 
PMMA resist (950PMMA C2 Resist from MICROCHEM) was spin coated at 
5000 rpm and baked at 180°C for 3 min. As prepared sample was patterned using e- 
beam lithography with dosage 200 |aC/cm^ and beam intensity 8. The exposed 
samples were developed using MIBK: IPA solution for 90 s and IPA for 30 s. The 
developed samples were imaged using scanning electron microscope (SEM by Tescan 
Vega), confirming that the periodic holes were exactly fitting the design. After that an 
Ag layer of 20 nm thickness was deposited directly on top of the sample using thermal 
evaporation of Ag pellets (from Kurt J. Lesker Company). As prepared samples were 
again imaged using SEM to see the uniformity of Ag layer on the gratings. Variable 
Angle Spectroscopic Ellipsometry (J. A. WooUam Co., Inc, V-VASE) was used to 
obtain the thicknesses and optical constants of Ag and AI2O3 thin films. The 
reflectivity spectra as a function of incident angle and excitation wavelengths were 
acquired using a J.A. WooUam Co. Inc V-VASE instrument with an angular 
resolution of 1° and wavelength spectroscopic resolution of 1 nm. The transmission 
measurements were performed using both J.A. Woollam Co. Inc V-VASE instrument 
and UV-Vis Spectrophotometer (Cary Series). 

Numerical simulations. The decay rate enhancement was numerically computed 
with Matlab. The outcoupling property of diffraction grating is introduced in the 
Fresnel reflection coefficients calculation, where the wavevector of the incident beam 
is taken as = nokosinO ± m(27r/AJ ± n{2n/Ay) with A^ = A^ = A for grating with 
square lattice symmetry. Here we assumed that + m and + n diffraction orders are 
responsible for modes excitation. However, commercially available lumerical FDTD 
software was used to simulate the electric field distribution. In the numerical 
simulation, the periodic boundary condition with smallest spatial grid size of 1 nm is 
used for the iteration to maintain the accuracy and stability of the FDTD simulation. 
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